Electroweak stars: how nature may capitalize on the standard model's ultimate fuel 
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We study the possible existence of an electroweak star - a compact stellar-mass object whose 
central core temperature is higher than the electroweak symmetry restoration temperature. We 
found a solution to the Tolman-Oppenheimer-Volkoff equations describing such an object. The 
parameters of such a star are not substantially different from a neutron star - its mass is around 1.3 
Solar masses while its radius is around 8km. What is different is the existence of a small electroweak 
core. The source of energy in the core that can at least temporarily balance gravity are standard- 
model non-perturbative baryon number (B) and lepton number (L) violating processes that allow 
the chemical potential of B + L to relax to zero. The energy released at the core is enormous, but 
gravitational redshift and the enhanced neutrino interaction cross section at these energies make the 
energy release rate moderate at the surface of the star. The lifetime of this new quasi-equilibrium 
can be more than ten million years. This is long enough to represent a new stage in the evolution 
of a star if stellar evolution can take it there. 

PACS numbers: 04.50.Kd 



CD 



> 
o 

o 

On 
O 



X 



INTRODUCTION 

The last stage of the evolution of a star massive enough 
to go through the supernova explosion (but not massive 
enough to become a black hole) is believed to be a neu- 
tron star. A neutron star has no active energy source 
which balances the force of gravity. Instead, it is sup- 
ported by degeneracy pressure, which arises due to the 
Pauli exclusion principle. However, this pressure can sup- 
port only neutron stars lighter than about 2.1Mq, which 
is known as the Tolman-Oppcnhcinicr-Volkoff limit. Neu- 
tron stars heavier than this limit eventually become black 
holes. 

Recent studies point out the existence of a state be- 
tween the neutron star and black hole. It is called a 
quark star. This state owes its existence to the QCD 
phase transition in which the original nuclear matter be- 
comes quark matter. This process can release 10^'^ erg in 
about 10"'^ to 10~^s [l|,01 in the form of neutrino bursts. 
While this a huge amount of energy, it is released in too 
short a time to provide the pressure to prevent gravita- 
tional collapse. After this energy release, a star contain- 
ing effectively only three quark flavors (u,d,s) can exist 
in a stable equilibrium where the pressure is provided by 
the Pauli exclusion principle. However, at higher densi- 
ties where four or more quark flavors are present, quark 
matter cannot avoid gravitational collapse. 

Since in the gravitational collapse matter gets com- 
pressed to ever increasing densities/temperatures, it is 
natural to explore what could happen at the electroweak 
phase transition, the next (and last) within the standard 
model. In the standard model, both baryon and lep- 
ton number U{1) global symmetries are accidental, con- 
served perturbatively, but violated by non-perturbative 



processes such as those mediated by instantons. At tem- 
peratures well below the electroweak symmetry-breaking 
scale, T <, lOOGeV, these non-perturbative processes are 
suppressed by the extremely small factor e"^'^/" (giving, 
for example, a proton life-time of lO^^-'^yrs Q). (Note, 
unless otherwise indicated we shall employ natural units 
where h = c = = i-) However, above this scale, 
baryon and lepton number violating processes are ex- 
pected to be essentially unsupprcssed (although the com- 
bination B — L remains conserved) . Quarks can then be 
effectively converted into leptons. In this process, which 
we call electroweak burning, huge amounts of energy can 
be released. 

Let us assume that in the core the temperature is 
above, or at least very near, the weak scale, and the mat- 
ter density is also comparable p ^ (lOOGel^)''. At these 
densities, the matter is opaque even to neutrinos and the 
energy released at the center cannot stream freely out 
of the star. Nevertheless, the energy can be carried out, 
mostly by neutrinos and anti-neutrinos. These can also 
carry out any excess anti-lepton number generated in the 
electroweak burning, and thus prevent the lepton num- 
ber chemical potential from halting the consumption of 
the baryon number. This mechanism can likely provide 
a stable energy source which can counteract gravity for 
a while. On a side track, we remark that studying the 
electroweak star physics, in particular the consequences 
of a very powerful energy source that was so far neglected, 
may provide important new clues about supernova explo- 
sions, but we reserve further discussion on this matter to 
a future paper (in preparation). 

In this paper, we study the inner structure of such a 
star and calculate its lifetime. We find that this new 
phase can last up to lO^^s, which is long enough to give 
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a new name to this type of stars - electroweak stars. In 
order to find out whether electroweak stars are very com- 
mon or very exotic objects, much more detailed analysis 
is required. Future work will consider questions of stabil- 
ity, of evolution - if and when the electroweak star arises 
in the late-phases of evolution of ordinary stars - of the 
structure of the outer cooler layers, and of observational 
signatures of these objects. 



DENSITY, PRESSURE, AND PARTICLE 
NUMBER DENSITY INSIDE THE 
ELECTROWEAK STAR 

We separate the core of the electroweak star into three 
regions, as shown in Fig. [TJ The central region is hotter 
and denser than the electroweak symmetry restoration 
temperature (T ^ lOOGeV). This region is very dense 
[p > lO^GeV^), but small (several cm). We find that 
the total mass inside this region is ~ 5 x 10~^Mq. Since 
the lepton and baryon number arc not conserved in this 
region, but B — L is, baryons are freely transformed to 
anti-leptons so long as it is thermodynamically favorable 
to do so. A S'f7(2)-preserving instanton interaction can 
convert 9 quarks into 3 anti-leptons for example: 

udd CSS tbb — >■ l/e , l'/^ , i^r (1) 

and all B — L and electromagnetic charge conserving 
variations involving quarks, anti-quarks, neutrinos, an- 
tineutrinos and charged leptons. At these densities and 
temperatures each particle carries about lOOGeV of en- 
ergy, so this process can release about SOOGeV per neu- 
trino. The typical energy of the produced neutrino is not 
the only relevant quantity here, we also need the rate at 
which the energy is released. The sphaleron baryogene- 
sis/baryodestruction rate (above the electroweak symme- 
try restoration temperature) is proportional to T"*. How- 
ever, we will show later that the energy release rate from 
the core is limited by the quark supply rate (dictated 
by the gravitational collapse) rather than the maximal 
electroweak burning rate. 

The energy in the form of neutrinos flows out of the 
central core, and eventually out of the star. However, 
in the central region, the mean free paths of all particles 
are short compared to the size of the star. The energy 
must therefore diffuse out of the central core. Similarly, 
the non-perturbative B-violating processes, while reduc- 
ing B, reduce L. Initially B — L > in stars, because 
neutrinos escape, and thus there are more neutrons plus 
protons than there are electrons. Thus, the B-violating 
processes would cause an anti-lepton number to build up. 
This must cither stop the burning, or be carried out of 
the core by a flux of anti-leptons. Outside the central 
core, baryon and lepton number arc conserved. The den- 
sity falls with increasing radius, and particle mean free 



paths increase, however, they are still too short for parti- 
cles to stream freely out of the star. Eventually, however 
we reach the star's neutrino release shell (i.e. neutri- 
nosphere), where the neutrino (and anti- neutrino) mean 
free path is long enough for these particles to flow freely. 
[This is well within the radius at which the same is true 
of photons.] The neutrinos and anti-neutrinos carry off 
not just most of the energy produced in the core burning 
region (the rest is carried off by photons), but also the 
lepton number, thus permitting the burning to continue. 
Since the fuel is predominantly baryons not anti-baryons, 
resulting in the production of an excess of anti-leptons 
over leptons, there will be a greater flux of anti-neutrinos 
than neutrinos from the star. Moreover, the electroweak 
processes (unlike inverse beta decays in normal stars) cre- 
ates equal numbers of anti-neutrinos of all three genera- 
tions (see Eq. [Ij . Thus the excess is likely to be equal in 
all generations. Such a signal, taking into account neu- 
trino oscillations, might be an indication that the elec- 
troweak processes are taking place at the core of the star. 
However, the excess is likely to be very small, since (as in 
a supernova), the emission is expected to be dominated 
by thermal emission from the ncutrinosphcrc (which car- 
ries no excess), not by direct anti-neutrino production 
in the core. More likely, a detailed analysis of the elec- 
troweak star will show that the neutrinospheres for the 
three generations of neutrinos and anti-neutrinos are sep- 
arated due to the differing masses and hence abundances 
and momentum distributions as a function of radius of 
the quarks and charged leptons of the associated gener- 
ations. Therefore the fluxes and spectra of the different 
neutrino and anti-ncutrino species should be predictably 
(and hopefully observably) different. 

The size of a quark star is about 10km, while its mass 
is of the order of Mq . A small electroweak core cannot 
significantly change the region outside the core. We use 
the so-called "bag model" [5[ as an approximate solution 
of the region outside the core. The pressure and energy 
density are Q 

P = J2p^-Bn (2) 

i 

e = + Bn- 

% 

Here, P is the total pressure, pi is the partial pressure 
of particle i, B^ is the bag energy density, e is the to- 
tal energy density, and is the energy density of par- 
ticle i. The subscript i runs over all types of particles 
present in the star (however here for simplicity we include 
only fermions, i.e. quarks, neutrinos and electrons). We 
choose the specific value B]^^ = 145MeV. The condition 
for the electric charge neutrality is 

qirii = 0. (3) 
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SOLUTION TO THE 
TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATIONS DESCRIBING AN ELECTROWEAK 

STAR 



FIG. 1: The structure of an electroweak star. We separate the 
core of the star into three regions: I. central core where the 
electroweak symmetry is restored, II. region where the high 
energy neutrinos are trapped, III. region from which high en- 
ergy neutrinos can escape. Neutrinos will be emitted between 
the regions II. and III. Photons are emitted at the surface of 
the star. 



Here, qi is the charge of the particle i. The pressure, 
energy density, and number density of particles can be 
well approximated from Fermi-Dirac and Bose-Einstein 
gas distribution: 



poo 

^/ dE{E'^mjy/^EME). 



(4) 
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Here, gi is the number of degrees of freedom of particle 
i (e.g. 5 = 2 for leptons and g ~ & for quarks), rrii is 
the mass of the particle, fjii is the chemical potential of 
the particle and T is the temperature. The distribution 
functions are 



h{E) 



for bosons and fermions. 
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1 ± e(^~^')/^ 



(5) 



In this section we setup and numerically solve the sys- 
tem of equations that describe the structure of an elec- 
troweak star. 

Outside the central region I, baryon and lepton num- 
bers are conserved. We assume for simplicity that the 
chemical potentials of the left- and right-handed fermions 
are the same, and that all fermions of the same charge 
(e.g. u, c and t quarks) have the same chemical poten- 
tial. An antiparticle has the chemical potential of the 
opposite sign. The interaction d ^ u + e + gives one 
relationship 



f^d= ^J.u + Me - Ml.- 



(6) 



Eq. Ilj implies another relation inside the electroweak 
core: 



(7) 



We assume fi,^ = outside the electroweak core, so that 
neutrinos or anti-neutrinos can propagate away. We 
set the boundary conditions at the center that pres- 
sure and energy density arc (lOOGcV)^, and then cal- 
culate the structure of the star according to the Tolman- 
Oppcnheimcr-Volkoff equations 



dP 
dr 

dM 
dr 

gtt 



{e + P){M + AirPr^) 
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2M(r) 



(8) 

(9) 
(10) 

(11) 



^^surface is radius of the star, Mgtar is the mass of the star, 
while Mp is the Planck mass. Since we are dealing with 
very dense objects, general relativistic effects described 
by the metric components gn and g,.,. will be very impor- 
tant. 

If we supply an equation of state e(P), the Tolman- 
Oppcnheimer-Volkoff equations could be solved directly. 
However, we do not want to assume an equation of state, 
and instead we try to solve the system of coupled equa- 
tions we already set-up. 

Because of the properties of Fermi-Dirac statistics, 
fermions can have large energies at T = 0. In that case 
the energy of neutrinos is not proportional to the temper- 
ature, instead it is proportional to the chemical potential 
as can be seen from Eq. ([5|). We will use this fact 
to eliminate T from the equations and simplify our cal- 
culations. The exact treatment would require keeping 



4 



T as a variable (beside the energy density, pressure and 
number density), setting the boundary condition at the 
center T = lOOGeV, and finding the temperature profile 
T(r). Setting T = will not significantly change P{r), 
e(r) and M{r), but may change n(r). This change in n(r) 
will seemingly aff'ect the neutrino mean free path, but one 
can explicitly verify that this does not happen since the 
mean free path in Eq. ([T3|) will effectively depend only 
on the energy density = niEi. Thus, the structure of 
the star will not change much if we set T = 0. [Note 
that the properties of Bosc-Einstcin statistics are differ- 
ent and one cannot set T ~ when including photons 
and weak gauge bosons into consideration.] 

Thus, we have four independent equations (jS]), ([6|). 
(O; and ([8|), with four independent variables ^u, ^J■d^ 
jjLe and ^u- The chemical potentials /id, /ie and /i^ 
then uniquely specify P{r), £{r), ni{r) and M{r) through 
equations ©, (U), ©, (H]) and dH). 

The numerical results are plotted in Fig. [2] The total 
size of the star is 8.2km (radius where pressure and en- 
ergy density drop to zero) and its mass is I.SMq (mass 
within that radius). The size of the central electroweak 
core is several cm (this is the distance at which the den- 
sity falls below (lOOGcV)^. Within this central region, 
the pressure and energy density gradients are very small, 
so pressure and energy density are practically constant 
there. We can also find the dependence e(P) and find 
that the equation of state does not change much across 
the star and is approximately e = 3P. In the absence of 
the active energy source, such an equation of state would 
signal instability. However, electroweak stars have a very 
powerful energy source which can balance gravity. Also, 
this equation of state will tip toward a non-relativistic 
equation of state once weak gauge bosons (and interac- 
tions among the particles) are included into calculation, 
since they are not ultra-relativistic at these energies. 

For comparison, in Fig. [3] we plot the energy density 
profiles of an electroweak star, a neutron star (a star 
made from neutrons only) and a polytropic star (with 
the polytropic index 1.5), where, to emphasize the differ- 
ence (and similarity) in the structures, we set the central 
density for all three to 3 x 10^°MeF^. (Of course, no 
stable neutron star can have this high a central density.) 
The neutron and electroweak stars are almost indistin- 
guishable at small radii because particle masses arc not 
important in that region, but they become distinguish- 
able at large radii. Furthermore, we are currently treat- 
ing the quarks in the electroweak star as a degenerate gas 
of non-interacting particles. Adding in the interactions, 
which is likely to improve the stability of the electroweak 
star, is the subject of a future work. 

The non-perturbative electroweak interactions provide 
the source of energy that keeps the star from collapsing. 
As the baryons are burned, gravitational collapse brings 
more material to the central part of the core feeding the 
fire. We assume that eventually a quasi-equilibrium be- 
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FIG. 2: The pressure, energy density and particle number 
density change with the radius of the star. The total radius 
of the star is about 8.2km (radius where pressure and energy 
density drop to zero), while the mass is about 1.3Mq (mass 
within that radius). 
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FIG. 3: Comparison between the energy density profiles of 
an electroweak star, neutron star (a star made from neu- 
trons only) and polytropic star (with the polytropic index 
1.5). For all three plots the central density is 3 x lO^^MeV*. 
The neutron and electroweak stars are almost indistinguish- 
able at small radii because particle masses are not important 
in that region, but they become distinguishable at large radii. 

twccn the burning and the refueling is achieved. If so, the 
star enters a quasi-stable state that can last for a while. 
The question of dynamical stability will be addressed in 
future work. 



LOCATION OF THE NEUTRINOSPHERE 

Neutrinos near the central region have very high ener- 
gies (but below the weak scale) and the local density is 
extremely high. Their mean free path in this region is 
thus very short, and they cannot stream freely from the 
star. However, as the local density falls with radius, and 
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the neutrinos' average energy decreases, their mean free 
path increases. We define the neutrino escape radius (as 
a function of neutrino energy) as the distance from the 
center at which a neutrino with that energy has a mean 
free path greater than the thickness of the overlying mat- 
ter and can therefore escape the star. The neutrinosphere 
is the set of such spherical shells for neutrinos of energies 
characteristic of the local temperature. 

Now we study the neutrino mean free path in dense 
media as presented in We assume that the cross 
section of neutrino scattering is 0] 



(12) 



Here, Gp 1.166 x 10~'^GcV^^ is Fermi constant, s is 
the center of mass energy, Ei is the energy of the particle 
i, while El, is the neutrino energy. The mean free path is 



1 = E'-^^ 



(13) 



Unlike any other active star, particles propagating 
through the electroweak stars suffer large gravitational 
redshift. Thus E^, is not a constant, but changes as 



EAr) 



E,, 



(14) 



The redshift ratio is shown in FiglH Obviously, the red- 
shift effect near the center is much larger than at the 
surface. Since inside the star e ~ 3P, the redshift can be 
estimated from Eq. ((TO)) . It is 



E,. 



Pjr) 
Pirn 



1/4 



(15) 



The pressure is about (lOOGeV)^ near the center, and 
about (lOOMeV)'' near the surface (before the bag energy 
changes the pressure a lot). Therefore, a particle retains 
only a fraction 10 ~^ of its original energy. Though in- 
stanton processes release huge amounts of energy at the 
center, the star can emit much more moderate amounts. 
The result is shown in Fig. \5\ A particle with the origi- 
nal energy of lOOGeV near the center indeed carries only 
165MeV as it leaves the surface due to the redshift. This 
energy will be further reduced due to opacity of the high 
density medium that neutrinos are traveling through. 

We now use this energy corrected for the redshift to cal- 
culate the mean free path on the surface. From Eq. (fT2|) 
and Eq. ([TB]) the mean free path on the surface is 



1 
A 



surface ) 



(16) 



Eq. (Ull) yields A - W^McV"'^ - 10~'*km. At the core, 
Eq. p6p yields A ~ 10~^^m. A neutrino with this energy 
cannot leave the star freely, and must lose its energy while 
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FIG. 4: The redshift factor \/ gtt{r) vs. radius of the star. The 
redshift near the center is much larger than at the surface. A 
particle with the original energy of lOOGeV near the center 
carries away only 165MeV as it leaves the surface. Opacity of 
the dense medium will further reduce this energy. 



propagating from the center to the surface. This energy 
loss will in turn further change the neutrino mean free 
path. Properly calculated neutrino energy (plotted in 
Fig. [71) will imply that the neutrino mean free path at 
the neutrinosphere is A '--^ 0.1km. 

With the redshift taken into account, the neutrinos 
have roughly the same energy as the background they 
are propagating in, which is true in general for the other 
types of particles. If a particle has the energy density e. 



its energy can be estimated as ~ e 



1/4 



P^^^. Comparing 



this with Eq. ([TSl) , the energy of a neutrino decreases with 
the distance the same way as the background. Thus, 
the neutrinos coming from the center and arriving to a 
certain point inside the star have roughly the same energy 
as the neutrinos that were emitted by thermal processes 
at that point. 



TRANSPORT OF ENERGY THROUGH THE 
STAR, ENERGY RELEASE RATE, AND THE 
LIFETIME OF AN ELECTROWEAK STAR 

We would like now to estimate the (neutrino) lumi- 
nosity of the electroweak star. A naive luminosity (i.e. 
energy release rate) can be obtained from the black-body 
radiation at the stars core 



Ci, < Cbb ^ (Tbhei/Tri? 



(17) 



This places an upper bound of about 3 x W^^MeV^ = 
7 X lO^^erg/s. At this rate it would take less than a 
second to release IMq. Fortunately, Eq. ([T7| is a severe 
over-estimate of the luminosity. It does not take general 
relativity into account, nor does it allow for the fact that 
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FIG. 5; Neutrino energy as a function of the distance from 
the center where only effects of gravitational redshift are taken 
into account (no energy loss due to interactions) . The energy 
on the surface is about 1000 times smaller than the original 
energy at the center. 

the luminosity depends not just on the temperature and 
density at the neutrino-sphere but also on their gradients, 
since those determine the net outward flux of energy. 

First we will estimate the energy release rate when we 
include redshift effects in calculating the ratio dE/dt. An 
upper bound on the energy release rate can be obtained 
from the free fall time of the incoming quark shell into the 
electroweak-burning core. This determines the maximum 
rate at which the electroweak core burning can be fed. 
The free fall acceleration over a short distance given by 
the mean free path A is roughly given by a = A/ew/'^ewj 
where row refers to a radius one mean free path (A) away 
from the electroweak core (relativistic corrections will not 
change the result significantly). Therefore, the energy 
release rate can be no more than 

fdE\ 2 / 2AAf(rew) . ^ . ^ ^ 

This is just lO^^MeV'^, which is already much lower 
than the earlier upper limits obtained from Eq. ([T7| . 
We compare this rate with the electroweak baryogene- 
sis/baryodestruction interaction rate, which is above the 
electroweak phase transition proportional to T'^ [l^l 

^^A.rtx20atT' (19) 

This rate is about lO'^^MeV^, which is much larger than 
the rate in (fT8)) . Therefore we can reasonably assume 
that as quarks reach the electroweak core, they are con- 
verted into neutrinos instantaneously. Otherwise, the in- 
falling matter would pile up and form a black hole. This 
result also implies that the star's structure and energy 



released at the surface do not crucially depend on pre- 
cise details of the electroweak processed at the core, i.e. 
the star is capable of burning much more quarks than it 
actually does. 

We now calculate the neutrino energy release rate at 
infinity more accurately taking into account the relativis- 
tic transport of energy through the star. Gravitational 
redshift and time delay along with the enhanced neutrino 
interaction cross section at these energies will affect the 
energy transport through the star and make the energy 
release rate moderate at the surface of the star. Imagine 
two surfaces inside the star whose radii are one A apart. 
The outward flux of neutrinos will be partially canceled 
by the backward flux due to neutrino collisions with other 
particles, and we are interested only in the net outward 
flux. The relativistic transport of neutrino energy can 
then be described by 

,_^,,diS{r)gu{r)) ^ ^ ^^0) 
\/-grrir)dr 

Here S{r) = dE{r)/{dtdA) is the energy flux at the ra- 
dius r, A is the area, while L = dE/dt is the energy 
release rate at inflnity (which is about the same order as 
the intensity at the neutrinosphere) . The factor 5tt(r) 
describes both the redshift and time delay. The factor 
drrif) takes into account that the space-time inside the 
star is not flat. The energy flux S{r) can be found from 
the energy density ei, given by (j4]) knowing the energy 
(i.e. chemical potential) of neutrinos E^,. The metric 
coefficients are given by Eq. ([TOl) and (fTTj). It is then 
straightforward to calculate L, given all the other pa- 
rameters. 

The neutrino mean free path is not a constant through- 
out the star, but changes from point to point. We there- 
fore set a practical definition of the location of the neu- 
trinosphere as i?nsph with 

where A(r) is given by Eq. (|16p . This definition requires 
that a neutrino, at a given radius r and a given energy E^, 
interacts only once with other particles before it leaves 
the star. 

With these definitions, we can calculate intensity L as 
the function of the radius of the neutrinosphere -Rnsph- 
From the condition (|2T|) . we find the neutrino energy i?^, 
which in turn gives the required energy flux S{r) for a 
given i?nsph. TheresultisplottedinFig.ini Larger radius 
-Rnsph implies lower energy density, which means that 
higher energy neutrinos can escape. The energy release 
rate therefore increases with the radius of the shell and 
it is maximal at the surface of the star. However, there 
it exceeds the limit from the quark free fall. This implies 
that the neutrino release shell must be inside the star. 
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FIG. 6: The energy release rate vs. neutrino release radius: 
The energy release rate increases with the radius of the neu- 
trino release shell (i.e. neutrinosphere) . The maximum en- 
ergy release rate is on the surface, but there it exceeds the 
limit from the quark shell free fall into the electroweak core 
(which minimizes the denominator in (|17|l 'l. This implies that 
the neutrino release shell must be inside the star. 

At the neutrinosphere, i?nsph, the following must be 
true, AT:r'^S{r = Rnsph)gttir = i?nsph) = L, since there 
is no backward flux of neutrinos (all the neutrinos that 
reach i?nspii leave the star freely). From there wc can find 
the neutrino energy, E^, at the central electroweak core 
needed to support a certain intensity at the surface of 
the star. Fig. [7] shows the results. Higher energy release 
rate needs the source of higher energy. If the neutrino 
release radius is 8.1km, then the original neutrino energy 
must be about SOOGeV, which is equal to the energy per 
neutrino that typical electroweak processes can release at 
the center. We therefore conclude that the radius of the 
neutrino release shell can be at most 8.1km. If the neu- 
trino escape radius is about this size, the energy release 
rate is about lO^^'McV^. Integrating this, we find that it 
takes about lO^^sec (10 million years) to release energy 
equivalent of IMq. This is the minimal life-time of the 
electroweak star, provided that the electroweak burning 
does not stop before all the quark fuel is spent. 

In our calculations we neglected the fact that some 
fraction of energy is carried away by photons, since they 
have shorter mean free path. Further, though the net 
lepton number is conserved in this region, the number of 
neutrinos and antineutrinos are not conserved. We also 
ignored the effects of energy transport due to convection. 
Finally, neutrino energy was calculated from the chemi- 
cal potential change, instead of the temperature change. 
However, we assume that these effects will not change 
the order of magnitude estimate. 

On a completely different side, it is not difficult to 
imagine that the electroweak physics could be responsi- 
ble for some of the Type II supernova explosions. 



FIG. 7: The neutrino energy at the surface of the central 
electroweak core vs. the radius of the neutrino release shell 
(i.e. neutrinosphere) i?nsph. From Fig. [6] we see that the 
energy release rate increases with the -Rnsph- Higher energy 
release rate needs the source of higher energy. The energy 
of neutrinos produced at the core therefore increases with 
the radius of release shell. However, the energy is already 
larger than 500GeV at 8.2km. The star cannot support this 
amount of energy implying that its radius must be smaller 
than that. In fact, SOOGeV neutrinos produced by baryon 
number violating processes at the core imply that the radius 
of the neutrino release shell is 8.1km. 



Many decades have elapsed since modeling of the core 
collapse supernova mechanism began in detail. However, 
the exact mechanism remains elusive so far. It appears 
that the explosion of supernova is a multiscale and mul- 
tiphysics event, where gravity, neutrino transport, fluid 
instabilities, rotation, magnetic fields etc are all very im- 
portant ingredients. The existing models failed to pro- 
duce the supernova explosion. It will be of high impor- 
tance to check if the electroweak physics can help super- 
novae explode by providing a very powerful energy source 
that was neglected so far in the studies of this problem. 
The expected final outcome of the stability analysis of 
the electroweak stars sequence will perhaps yield a sta- 
bility curve. The possible outcomes are that the whole 
electroweak sequence is stable or unstable. However, a 
generic outcome with different regions in parameter space 
corresponding to stable and unstable stars depending on 
the initial conditions is more likely. In unstable regions, 
the electroweak stars with a certain initial mass and cen- 
tral density could be unstable to explosions. These ex- 
plosions could be in the form of gamma ray bursts or 
perhaps supernova explosions. To find out the answer, a 
very detailed modeling of the electroweak star is needed 
taking all the necessary ingredients like gravity (with full 
GR effects), neutrino transport, fluid instabilities etc into 
account. 
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CONCLUSIONS 

We studied the possibility of the existence of a new 
phase in the stelar evolution. After all of the thermonu- 
clear fuel is spent, and possibly after the supernova ex- 
plosion, but before the remaining mass crosses its own 
Schwarzschild radius, the temperature of the central core 
may surpass the electroweak symmetry restoration tem- 
perature. At this point, non-perturbativc baryon-number 
violating interactions in the ordinary SU{2) x U{1) elec- 
troweak standard model become unsuppressed. The con- 
sequent relaxation of the stellar baryon number chemical 
potential, i.e. the conversion of baryons to anti-leptons, 
may provide a source of pressure which can balance grav- 
ity. We constructed a solution to TOV equation whose 
central pressure is non-singular (and thus not a black 
hole). We showed that the in- falling matter gets con- 
verted into neutrinos at the rate much faster than the 
free-fall rate, which indicates that the mater burns be- 
fore crossing its own Schwarzschild radius (assuming that 
the core is not within the Schwarzschild radius itself). 
Our analysis shows that lifetime of this new phase can 
be as long as 10 million years which we propose to call 
an electroweak star. We emphasize that the electroweak 
star relies only on the Standard Model physics for its 
existence. 

Electroweak stars would be an exciting addition to the 
diverse menagerie of astrophysical bodies that the uni- 
verse provides (see also recent proposals in [lH). Never- 
theless, considerable work remains to be done before we 
can claim with confidence that such objects will form in 
the natural process of stellar evolution, or that they will 
indeed burn steadily for an extended period. Similarly, 
assessing the visibility of these fascinating new objects 
requires a careful modeling of their outer structure to 
determine the neutrino/ antincutrino and photon lu- 
minosity and spectrum. 

While we constructed a non-singular solution to TOV 
equations which describes the electroweak star, the cru- 
cial next step would be to verify that a full time de- 
pendent evolution will take the collapsing system to this 
configuration. This is not unlikely to happen since the 
global parameters of the electroweak star are not much 
different from those of the neutron star. The radius is 
8.2km and the mass is 1.3 Solar masses. What is differ- 
ent is the central core which is at or above the electroweak 
density. In gravitational collapse, the collapsing object 
is not crossing its own Schwarzschild radius at once as 
a whole. The simulations show that the central density 
is much higher than the rest of the star (since the outer 
layers push the inner layers stronger and stronger). It 
appears that it is not difficult to achieve central density 
of the order of the electroweak phase transition. For ex- 
ample, in [l3| . the authors find that densities achieved in 
gravitational collapse are much larger than what would 



be expected from the selfsimilar solutions (homologous 
collapse, PPN and other commonly used approximations 
fail in this regime). In particular, the Fig. 6 in [l3j shows 
that, with appropriate initial conditions (say those of a 
neutron star), the electroweak densities can be reached 
before the " gravothermal catastrophe", i.e. black hole 
formation. Once that stage is reached, electroweak burn- 
ing will get ignited and the collapse will be stopped for a 
while. This certainly does not imply that the whole star 
is at the electroweak density, in contrary, it is not much 
different from the ordinary neutron star, as the TOV so- 
lution demonstrates. 

One intriguing possibility is that the authors of 
and 15[ are correct that gravitational collapse of a stel- 
lar mass does not result in a black hole but in an object 
of very high density and temperature toward its center. 
Then one would not have to worry about the collaps- 
ing object crossing its own Schwarzschild radius before 
reaching electroweak densities. 

Finally, electroweak phase transition may be triggered 
inside a neutron star by a mechanism reminiscent of sono- 
luminescence, where non-linear waves (e.g. star-quakes) 
may concentrate huge energy within some small volume 

m 

We emphasize further that electroweak stars do not 
have to be formed exclusively in stellar matter collapse. 
Primordial black holes are formed in large density fluc- 
tuations in the early universe, with masses ranging from 
one Planck mass to a few solar masses, which is the mass 
within the horizon at the QCD phase transition (phase 
transitions usually produce fluctuations large enough to 
cause the collapse of the whole causally connected re- 
gion). Electroweak densities are reached very easily in 
such violent processes where usual homologous approx- 
imation clearly fails. In that case, however, the elec- 
troweak stars should be very long lived in order to be of 
astrophysical significance. One more caveat is of course 
that so far no primordial black holes have been observed. 

Studying the details of the the electroweak star 
physics, may give some clues for understanding of the 
supernova explosion mechanism and phenomenology on 
the side track. Also, the electroweak processes give an 
effective mechanism of erasing baryon number from a col- 
lapsing star before it becomes a black hole, which may 
be interesting from the point of view of the information 
loss paradox (see discussion in 17|). 

Seeing the gravitational collapse as an inverse Big Bang 
process, one might ask what happens when the next 
phase transition is hit, i.e the Grand Unifying Theory 
(GUT) phase transition. Our preliminary estimates in- 
dicate that the core of such an object should be micro- 
scopic. It is very unlikely that this stage will be reached 
before the object crosses its own Schwarzschild radius. 

DS acknowledges the financial support from NSF, 
grant number PHY-0914893. GDS and AL thanks CERN 
for its hospitality. GDS is supported by a grant from 



9 



the US DOE to the particle astrophysics theory group at 
CWRU. 



[1] A. Drago, A. Lavagno and I. Parenti, Astrophys. J. 659, 

1519 (2007) arXiv:astro- ph/051265 2| . 
[2] A. Drago, G. Pagliara, G. Pagliaroli, F. L. Villante and 

F. Vissani, arXiv:0809.0518 [astro-ph]. 
[3] F. C. Adams and G. Laughlin, Rev. Mod. Phys. 69, 337 

(1997) arXiv:astro-ph/9701131 . 
[4] N. S. Manton, Phys. Rev. D 28, 2019 (1983). 
[5] A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn and 

V. F. Weisskopf, Phys. Rev. D 9, 3471 (1974). 
[6] C. Kettner, F. Weber, M. K. Weigel and N. K. Glenden- 

ning, Phys. Rev. D 51, 1440 (1995). 
[7] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 

374 (1939). 

[8] C. Shen, U. Lombardo, N. Van Giai and W. Zu o, Phys. 

Rev. C 68, 055802 (2003) [arXiv:nu ci-th/0307101 . 
[9] W. J. Marciano and Z. Parsa, J. Phys. G 29, 2629 (2003) 



[10] 
[11] 



[12] 
[13] 
[14] 
[15] 
[16] 
[17] 



[arXiv:hep-ph/04 03168\ 
G. Moore, NucL 



[arXiv:hep-ph/9810313 



Phys. B 568, 367 (2000) 



D. Spolyar, K. Freese and P. Gondolo, Phys. Rev. 
Lett. 100, 051101 (2008) |arXiv:0705.052T1 [astro-ph]]. 
K. Freese, C. Hie, D. Spolyar, M. Valluri and P. Bo- 
denhei mer, arXiv: 1002.2233 [astro-ph. CO]. K. Freese, 
larXiv:0 812.4005 [astro-ph]. D. C. Dai and D. Stojkovic, 
JHEP 0908, 052 (2009) tarXiv:0902.3662. [hep-ph]]. 
S. Horiuchi, J. F. Beacom and E. Dwek, Phys. Rev. D 
79, 083013 (2009) ;arXiv:0812.3157 [astro-ph]]. 
S. Balberg, S. L. Shapiro and S. Inagaki, Astrophys. J. 
568, 475 (2002) arXiv:astro-ph/011056i; . 

C. Barcelo, S. Liberati, S. Sonego and M. Visser, Phys. 
Rev. D 77, 044032 (2008) [arXiv: 0712. 11301 [gr-gc]]. 

M. Visser, C. Barcelo, S. Liberati and S. Sonego, 
larXiv:0902.0346l [gr-qc] . 

W. Simmons, J. Learned, S. Pakvasa and X. Tata, Phys. 
Lett. B 427, 109 (1998) arXiv:astro-ph/9602143 . 

D. Stojkovic, F. C. Adams and G. D. Starkman, Int. J. 
Mod. Phys. D 14, 2293 (2005) 



arXiv:gr-qc/0604072 



